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Abstract 

Background: The effect of breed on Porcine Reproductive and Respiratory Syndrome Viremia (PRRSV) was tested 
using data collected in 17 Italian commercial pig farms and 1096 genotypes obtained by the PorcineSNP60 
BeadChip. A binomial logistic model was used to investigate the relationship between breed-clusters and PRRSV 
susceptibility. Breed-clusters were defined using the matrix of genomic kinship between all pairs of piglets. 

Results: Only the contemporary group effect, defined as all piglets reared in the same herd, in the same year and 
whose samples were collected in the same season, was significant. Sex, age and breed-cluster showed no 
statistically significant effect on PRRS viremia, although the Landrace and Cross breed-clusters showed the lowest 
Odds-Ratio 

Conclusions: The model failed to detect a significant breed-cluster effect, highlighting the impact of environment 
and management on PRRS viremia incidence. Incomplete exposure over the observed period may have masked 
possible breed differences. 



Background 

Over the last decade the genetics of host susceptibility/ 
resistance to disease has been a major topic of research. 
For the swine industry PRRS represents one of the most 
economically important diseases worldwide, causing 
reproductive failure, abortions, stillbirths, interstitial 
pneumonia and decreased growth rate [1]. Previous stu- 
dies investigated the possible role which breed effect 
may have in determining resistant/susceptibility of pigs 
to PRRSV. One author [2] in an experimental "in vivo" 
study found difference between Duroc, Meishan and 
Hampshire while two authors [3,4] with an "in vivo" 
infection experiment using Large White I Landrace and 
Hamphshirel Duroc synthetic lines found the former 
being more resistant to the effects of the virus. Ait-Ali 
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et al[5] studied the innate immune response to PRRSV 
infection "in vitro", using used flow cytometry to analyse 
cells in bronchial alveolar lavage fluid (BALF) from five 
commercial pig lines. Macrophages from the Landrace 
line showed significantly reduced virus replication and 
poor growth of PRRSV than Large White, Pietrain and 
other two synthetic pig lines. Whilst "in vivo" and "in 
vitro" experiment can possibly help in dissecting PRRS 
pathogenesis among breeds, these studies were usually 
based on relatively small numbers of animals. Moreover, 
they do not consider environment and management 
which actually play an important role in the host resis- 
tant/susceptibility to PRRSV. Using data collected in 
commercial pig farms and genotypes obtained by the 
PorcineSNP60 BeadChip (Illumina, San Diego, CA), the 
objective of the present study was to cluster animals 
based on the average similarity among them and to test 
a possible breed-cluster effect on PRRS viremia. 
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Methods 

The data for the present study were extracted from an 
existing database belonging to the MISAGEN project 
[6], which included pedigree information, clinical symp- 
toms, and health related phenotypes collected from a 
commercial pig breeding population in the north of 
Italy. The original dataset included records for the PRRS 
viremia measured by PCR in sera of 2908 weaning pig- 
lets from four breeds, namely LargeWhite, Landrace, 
Duroc and Pietrain. PRRS viremia was defined as a bin- 
ary trait based on the results of the traditional PCR: 
negative samples were coded as 0, positive samples as 1. 
DNA samples from 1,096 piglets were genotyped for 
64,232 SNPs using the PorcineSNP60 BeadChip [7]. 
Negative and positive piglets numbered 766 and 440 
animals, respectively. The number of negative/positive 
piglets by breed is shown in table 1. Prior to statistical 
analyses genotype quality control was performed using 
the quality control function implemented in the GenA- 
BEL package for R statistical software [8]. The following 
filters were applied for exclusion of individual single 
nucleotide polymorphisms (SNPs): 

- call-rate (%) < 99 (if the SNP was available in less 
than 99% of all genotyped individuals) 

- minor allele frequency (MAF) in all individuals < 0.05 
Furthermore individuals <99% call rate (maximum 

percent of missing genotypes in an individual) were 
eliminated. A total of 14,967 SNPs (24.8%), from the 
available 60,123 SNPs, were excluded from the case-con- 
trol analysis if one of the filters indicated a violation of 
the quality. A total of 77 (0.063%) individuals were 
excluded: 33 individuals who had low call rates, 3 indivi- 
duals with too high Identity By State (IBS) and 41 indi- 
viduals with sex discrepancies. A contemporary group 
(Herd-Year-Season, HYS) was defined as all piglets 
reared in the same herd, in the same year and whose 
samples were collected in the same season. Sampling 
season was categorized as season 1 (January to April), 
season 2 (May to August) and season 3 (September to 
December). Totally, 46 HYS groups were defined whose 
average size was 9.1. Three breeds were present in 
twenty HYS (43 %), two breeds in sixteen HYS (35%), 
four breeds in three HYS (15 %) and only one breed in 
seven HYS (15 %). A General Linear Mixed Model with 

Table 1 Number of negative (0) and positive (1) piglets 
by breed (original dataset) 



Breed 


0 


1 


% Infected 


Duroc 


189 


137 


42% 


Landrace 


134 


67 


33% 


Large White 


309 


207 


40% 


Pietrain 


24 


29 


55% 


Total 


656 


441 


40% 



lbs: ■ = — "V - 



where L is the 



a logit link function and binomial distributions was used 
to investigate the relationship between breed and 
PRRSV susceptibility using the Glimmix procedure in 
SAS version 9.1 (SAS Institute Inc., Cary, NC). Breed- 
clusters were defined using the matrix of genomic kin- 
ship between all pairs of piglets computed as: 

L {gi,i-Pi){gi,j-Pi) 
PlU-Pl) 

number of loci, pi is the allelic frequency at 1-th locus 
and gi ; j is the genotype of j-th individual at the 1-th 
locus, coded as 0, 1=2, and 1, corresponding to the 
homozygous, heterozygous, and other type of homozy- 
gous genotype. Such a matrix was transformed to a dis- 
tance matrix and the first two principal components 
were computed by Classical Multidimensional Scaling. 
The fitted model was: Y^i = HYS t + G ; + B k + F t + b. 
agei+ eij k i (1), where Y is the binary trait, HYS is the 
fixed effect of the ith contemporary group, G is the 
fixed effect of jth gender class, B is the fixed effect of 
the kth breed-cluster class, F is the random effect of 1th 
piglet, age is the covariate for age at sampling, with b 
being the regression coefficient, and e is the random 
error. 

Results 

Breed-clusters as a result of "a posteriori" definition 
based on Kinship analysis are plotted in Figure 1. Each 
point represents a piglet. As expected, using SNP infor- 
mation, four main clusters can be identified. Those clus- 
ters matched with the expected breed, defined by the 
time of field data recording. However, an additional 
cluster (cluster 4) can be identified between cluster 1 



Cluster 1=Large White 
Cluster 2=Pietrain 
Cluster 3=Landrace 
Cluster 4=Cross 
Cluster 5=Duroc 
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Figure 1 Multidimensional Scaling Plot - Breed Clustering. 
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Table 2 Number of negative (0) and positive (1) piglets 
by breed-cluster (after SNP editing) 



Breed-Cluster 


0 


1 


% Infected 


1 (Large White) 


271 


182 


40% 


2 (Pietrain) 


26 


28 


52% 


3 (Landrace) 


121 


63 


34% 


4 (Cross) 


10 


10 


50% 


5 (Duroc) 


180 


128 


42% 


Total 


608 


412 


40% 



(Large White) and cluster 5 (Duroc ). Those animals are 
possibly the result of cross-breeding. Indeed, they belong 
to farms where both breeds (i.e, Large White and Duroc) 
were present and where artificial insemination was used. 
Our hypothesis is that some inconsistencies occurred 
during the insemination events, i.e. that a limited num- 
ber of boars was associated to the wrong insemination. 
According to the breed-cluster classification (Table 2), 
the observed proportion of affected piglets ranges from 
34% (Landrace - Cluster 3) to 52% (Pietrain - Cluster 
2). Results from Type III test of fixed effects are shown 
in table 3. Only the Herd-Year-Season effect (contem- 
porary group) is significant. Sex, age and breed-cluster 
show no statistically significant effect on PRRS viremia, 
although the Landrace and Cross clusters show the low- 
est Odds-Ratio (not shown). Results on breed-cluster 
effect are contradictory. Some studies reported the exis- 
tence of differences in susceptibility among breeds 
[9,10]. Nevertheless, most of those studies were based 
on in vitro or in vivo experiments, used commercial 
data from few large farms or analyzed different pheno- 
types. The present study spanned a 3-year period and 
used data collected from 17 farms. As Cooper et al. [11] 
observed, data under controlled experimental conditions 
do not necessarily support field reports. Physical and 
environmental factors may affect the immune system 
and hence determine or change the host response. This 
is exactly what our results suggest, because the contem- 
porary group effect tries to capture the combined effect 
of environment and management removing variation 
due to their conditions over time. Previous findings 
from MISAGEN data [12] do show the existence of 
genetic variation but no breed effect. In this case breed 
was not defined by Kinship analysis. Bishop and Wool- 
liams [13] showed how incomplete exposure to infection 

Table 3 Type III text of fixed effects of breed, age, HYS 
and sex 



Effect 


Num DF 


Den DF 


F Value 


Pr > F 


Breed-Cluster 


4 


317.6 


0.91 


0.4597 


age 


1 


651.7 


0.24 


0.6216 


hyseason 


45 


172.5 


1.9 


0.0018 


sex 


1 


967 


0 


0.9815 
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Figure 2 Observed PRRS viremia incidence over the test period. 



can reduce the power of datasets. This is the case of the 
present study, where records were collected in commer- 
cial farms spanning a 3-year period. As can be observed 
in Figure 2 incidence varies from a minimum of 27 % 
(September - December 2007) to a maximum of 80 % 
(September - December 2008). 

Conclusions 

The idea behind this study was to use SNP to correctly 
cluster animals based on the average similarity among 
them and to test a possible breed-cluster effect using a 
General Linear Mixed model. The model fails to detect 
a significant breed-cluster effect but highlighted the 
impact of environment and management on PRRS vire- 
mia incidence. Although we cannot formally exclude 
that incomplete exposure over the observed period may 
have masked possible breed differences, the genome 
wide analysis currently in progress could detect a signifi- 
cant genetic variability in host response to PRRSV in the 
same dataset (Botti, Biffani et al. in preparation). 
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